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Abstract: Applying the AdS/CFT correspondence, the expansion of the heavy-quark potential of the N - A 
supersymmetric Yang-Mills theory at large Nc is carried out to the sub-leading term in the large 't Hooft 
coupling at a nonzero temperature. The strong coupling corresponds to the semi-classical expansion of 
the string-sigma model, the gravity dual of the Wilson loop operator, with the sub-leading term expressed 
in terms of functional determinants of fluctuations. The contribution of these determinants are evaluated 
numerically. 
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1. Introduction 

A prominent implication of the AdS/CFT duality ^, |^ is the correspondence between the type IIB 
superstring theory formulated on AdSs x and N - A supersymmetric Yang-Mills theory (SYM) in four 
dimensions. In particular, the supergravity limit of the string theory corresponds to the leading behavior of 
SYM at large A^e and large 't Hooft coupling 

A^g\^Nc^^. (1.1) 

with L the AdS radius and a' the reciprocal of the string tension. This relation thereby has been applied 
fruitfully to explore the strongly coupled QGP created in RHIC in spite of its underlying dynamics, QCD, 
is different from N - A SYM. See for a review and the references therein. It is expected some of the 
properties of the latter is universal for all strongly interacting systems. The heavy quark potential addressed 
in this paper is one example of such an expedition. It is expected that some discoveries of N - A- SYM at 
large Nc and large A may shed new light on real QCD and other strongly interacting systems (For a recent 
comparison with the lattice simulation of the large Nc QCD, see [|[)). 

The heavy quark potential of QCD is an important quantity that probes the confinement mechanism in 
the hadronic phase and the meson melting in the plasma phase. It is extracted from the expectation of a 
Wilson loop operator 

WiQ^Pe-'^cd^i-^' (1.2) 

with the gauge potential and the symbol P enforcing the path ordering along the loop C. The thermal 
expectation value < W{C) > can be measured for QCD on a lattice and the heavy quark potential is defined 
with F-ansatz or U-ansatz. The F-ansatz accounts for the induced free energy of a heavy quark and a heavy 
antiquark separated by a distance r 

J7r ■r^ ti < trW(C2) > 

F(r,T) = -T\n r, (1.3) 

^ ' <trW(Ci)>2 

where Ci is a straight line running in the direction of the Euclidean time, closed by the periodicity of the 
Matsubara formulation, and C2 consists of two such lines separated by r and running in opposite directions. 
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The self-energy of each heavy quark has been subtracted via W{Ci). The U-ansatz accounts for the induced 
internal energy and is given by 

^/(r,^) = -^2A^|J. (i.4) 

We have F(r, 0) = U(r, 0). 

In the case of AT = 4 SYM, the AdS/CFT duality relates the Wilson loop expectation value to the path 
integral of the string-sigma action developed in of the worldsheet in the AdSs x bulk spanned by the 
loop on the boundary for T - and its generalization in the presence of a black hole for T > To 
the leading order of strong coupling, the path integral is given by its classical limit, which is the minimum 
area of the world sheet, and the heavy quark potential of F-ansatz extracted for the loop C2 and Ci reads [|[ 

[|oii[[ni] 

F(r,T) = -—--—mm[go(rT),0], (1.5) 
r4(i) r 



where go(rT) is a monotonic decreasing function and is given by ( |2.14| ) of the next section. We have goiO) = 

).7541 

7lT 



1 and goiroT) - with the screening length ro - . The leading order potential of the U-ansatz is 



obtained by substituting (1.5) into ( p,.4j ) 



The strong coupling expansion of the SYM Wilson loop corresponds to the semi-classical expansion of 
the string-sigma action and is parameterized by 

with gi(0) = 1. The coefficient k, corresponding to the subleading term at T = 0, has been calculated 
numerically in and analytically in [p3[]. We have k ^ -1.33460. Computing the function gi(rT) is the 
main scope of this work. We find that g\{rT) is a monotonically decreasing function, which reaches 0.92 
at ro where goiroT) = 0. Therefore the temperature correction to the subleading term of the heavy quark 
potential is numerically small. The screening radius to the subleading order, r'^, determined by 



VI 

reads 



go(roT)+ " ^ ^0 (1.7) 
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ro = ro . (1.8) 

VIr 



This work is the continuation of the previous calculation at zero temperature [|2|. The classical solution 
of the string-sigma model and the one loop effective action underlying k is briefly reviewed in the next 
section. There we also outline our strategy of computation, which is within the framework of [ p^ and [p^. 
The computation of the functional determinants involved by Gelfand-Yaglom method is described in 
the sections III and IV. The section V concludes the paper along with some discussions of the result and 
some open issues. 

2. The formulation 



Let us begin with a brief review of the classical limit that leads to the leading order potential (1.5). The 
string-sigma action in this limit reduces to the Nambu-Goto action 



with g the determinant of the induced metric on the string world sheet embedded in the target space, i.e. 

dX'' dX 
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where X'^ and Gf^y are the target space coordinates and the metric, and cr" with ff = 0, 1 parameterize the 
world sheet. The target space here is Schwarzschild-AdSs xS^, whose metric may be written as 



f 



(2.3) 



where / = 1 - §r with zi, - and dO.^ is the element of the solid angle of S^. The physical 3-brane resides 
on the AdS boundary z - 0. The string world sheets considered in this paper are all projected onto a point 
of in the classical limit. 

The Wilson loop of a static heavy quark, denoted by Ci, is a straight line winding up the Euclidean time 
< f < j0 periodically at the AdS boundary, where /5 = 1 /T. The corresponding world sheet in the AdS bulk 
can be parameterized by f and z with x constant and extends all the way to Schwarzschild horizon, z ^ Zh- 
The induced metric is given by 



ds'[Ci] = \(fdt^ + ^dz') 



with the scalar curvature in a diagonal metric 



R 



1 a — 
y[gdz 



We find 



(1) 



(2.4) 



(2.5) 



(2.6) 



Substituting the metric (2.4) into ( 2. 1 ), we find the self-energy of the heavy quark 

1 r~'' dz 

2na' Js Z^ 



(2.7) 



Notice that we have pulled the physical brane slightly off the boundary with the radial coordinate z = e, as a 
regularization of the divergence pertaining the lower limit of the integral ( 2.7). 

The free energy of a pair of a heavy quark and a heavy antiquark separated by a distance r can be 
extracted from the Wilson loop consisting of two anti-parallel lines each winding up the Euclidean time at 
the boundary. This Wilson loop will be denoted by C2 and the world sheet in the bulk can be parameterized 
by t and z with x\x^ - const. The function x^(z) is determined by substituting the induced metric 



(2.8) 



into the action (2.1 ) and minimizing it. We have 



^(4-z'^)(4-z'') 



(2.9) 



The maximum bulk extension of the world sheet, zo < Zh, is determined by the distance r between the two 
lines at the boundary and we find that 



zq Jo 



ed( 



V(l-^4)(,4_^4) 



(2.10) 



where s - z^. Substituting (2.9) into (2.8), we end up with the induced metric 



fdt^ + 



zt-z' 



dz^ 



(2.11) 
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and the scalar curvature 



1 - / 4/0 
R^-2 + 2- — U3 --^) = R^^^ 
1-/0 



(2.12) 



where /o = 1 - ^. The free energy of the heavy quark pair is given by 

^^fdz 



E[C2\ = TS'iCi] 



(2.13) 



where the same regularization is applied to the lower limit of the integral. 

The function goirT) of (1.5) that underlies the heavy quark potential is obtained by subtracting from 



(2.13) the self energy of each quark(antiquark), i.e. 

-T—goirT) = lim(£[C2] -2£[Ci]) = 

p4 / 1. j r £->o+ 



dz 



•Jzo ^ 



(2.14) 



and is divergence free. 

The one loop effective action, W is obtained by expanding the string-sigma action of |^ to the quadratic 
order of the fluctuating coordinates around the minimum area and carrying out the path integral [Q. We 
have 

det4(-,y«V„ + T3) 1 



W[Ci] 

for the static quark or antiquark and 

W[C2] = -ln[ 
for the quark pair, where 



.i„[ 



deti(-V2 + I + Rm)deti(-V^y 



det4(-/r«V„ + T3) 



det^ (- V2 + 4 + _ 26) det(- + 2 + 5) det^ (- V^) 



^= ^(7-l)(/-/o) 

4 f 



(2.15) 



(2.16) 



(2.17) 



The determinants in the denominators of (2.15) and (2.16) come from the fluctuations of three transverse 
coordinates of the Schwarzschild-AdS sector and five coordinates of with the Laplacian given by the 
metric (2.4) or ( ^.11 ). The determinants in the numerators come from the fermioninc fluctuations, where we 
have introduced the 2d gamma matrices, 70 = = T2, 71=7'= ti and 7071 - -irj, with ti , T2 and T3 the 
three Pauli matrices. In terms of the zweibein of the world sheet, e^, we have y^ = eiyj with y = 0, 1 and 
the CO variant derivative 

v„ = — + -[r/,r>K (2.18) 

dxa 8 



with (jjjj the spin connection corresponding to (2.4) or ( 2. 11 ). The power "4" comes from eight 2d Majorana 
fermions each of which contributes a power 1/2. The one loop correction to the heavy quark potential is then 



Fi(r, T) = T lim (WiCi] - 2W[Ci]). 



(2.19) 



The simplicity of the static gauge (Nambu-Goto action) adopted here is not cost free. The efi'ective 
action W[Ci] or W[C2] logarithmically divergent, the coeflicient of which is proportional to the volume part 
of the Euler character 



/ 



dtdz '\fgR 



(2.20) 
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of each world sheet with the same coefficient of propoitionaHty [y_6|. It follows from ( [2.4[ ), ( [2.6[ ), ( |2.11[ ) and 
( 2.12 ) that the 0{ -) term of the integral ( Z.20 ) for the parallel lines is exactly twice of that for the single fine. 



We have indeed that 



for the straight line and 



cf(Tyf^R=P^ -2 [-2(1 



3^)] 



S Zh 



(2.21) 



dz- 

z 



(2-/) 



1 - f 4/o 
[-2 + 2— ^(3- -f)] 

1-/0 



s 



(2.22) 



for the parallel Unes. Unlike the zero temperature case, however, the logarithmic divergence does not cancel 
in the combination of ( 2.19| ) in the limit e — > because of the 2nd term on RHS of (2.21). This may be 



the artifact of different topologies of C2 and two Ci 's in the presence of the back hole horizon. The integral 



(2.22) does not match twice of (2.21) as zq — > Zh, when the metrics (2.4) and (2.11) agree. Nevertheless, 



a finite one-loop correction to the heavy quark potential can still be extracted by the following argument, 
which is slightly handwaving 1 14|. Assuming that the one-loop formulation of [ [T^ in the Polyakov gauge 
can be generalized to the nonzero temperature case and the effective actions computed in this gauge, denoted 
by 'WlCi] and ^[Ci], are both finite and the one-loop correction reads 



Fi(r, T) = Ti^VlCi] - 2^[Ci]). 



(2.23) 



It was shown in [ |16| ] that ^[Ci] = at T = 0. Because T is the only dimensional quantity for ^[Ci], 
it must remain zero for all T . Next, we introduce a reference temperature To and the corresponding world 
sheet with the superscript "(0)", we have 



Fi{r, T) - Fiir, To) = T^[C2] - To^[Cf]. 



(2.24) 



Then, we identify this difference with that computed in the static gauge since the latter is free from the 
logarithmic divergence and is finite in the limit e — > 0. Since F(r, 0) has been obtained before [p^, all we 
need to calculate here is the difference Fi(r, T) - F\{r, 0). 



Parallel to the zero temperature case, we transform the metric ( 2.11 ) into conformal one 

ds\C2\ = e^Hdt^ + dcr^) 



(2.25) 



in order to remove the artificial singularity at z - Zo introduced by the double coverage of the world sheet C2 
on the coordinate patch (f , z), where cr is related to z via 



We have cr goes from -K to and then from to K as z moves from to zo and back to 0, where 

d^ 



K = \o-,„aA = ZqS- 



f 

Jo 



V(l-^'*)(,4_^4) 



(2.26) 



(2.27) 



Either zo or K serves a length scale of the problem. The first two terms of the expansion of z in the 
powers of (K + cr) reads 

Z^{K±a)[\-^{- + -){K±cTf] (2.28) 
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for ± cr <K zo- 



The conformal factor of (2.25) is given by 



/ 



(2.29) 



expressed in terms of cr by inverting ( [2.26| ). 

In a conformal metric,the scalar curvature R reads 



,01 



R 



-2e- 



2z^ da/; 
f do- 



where dLtj" is the spin connection 
To proceed, on writing 



we have 



' dcr z f dcr 



-V2 = e-2^Ao[C2] 
- + 2 + 5 = e-^^Ai [Ci] 
-V2+4 + 7;-25 = e"2%[C2] 



Ai[C2] 



+ e^'t'il + 6), 



and 



, d ldd)\ d . 



-ITl- iT2 — + T^e" 

OCT at 



(2.30) 

(2.31) 

(2.32) 
(2.33) 
(2.34) 
(2.35) 

(2.36) 
(2.37) 
(2.38) 

(2.39) 



As was shown in [p^, the scaling factors e ^'^,e and e=^2 do not contribute to the determinant ratio (2.16) 
and it becomes 



1V[C2] =-hl[- 



det4D^[C2] 



-] = -,„[ 



'detJAo[C2] detA2[C2] det^ Ai[C2] 
where we have used the relation 

det^Df [C2] = detD|[C2] = detA+[C2]detA-[C2] 

with 



det2A+[C2]det2A_[C2] 
detUo[C2] detA2[C2] detUi[C2] 



\, (2.40) 



(2.41) 



(2.42) 



Making a Fourier transformation of the time variable f, each functional determinant of (2.15) and (2.16) 
is factorized as an infinite product of its Fourier components with each Fourier component obtained by 
replacing the time derivative ^ in eqs. (2.36)-(2.39) by -ioj with the Matsubara frequency cj - 2nnT for 
the bosonic determinants and w - (2n + l)7rT for the fermionic ones. We have then 

00 

WlCi] = - 2 ^ [Zh[2nnT]Zf((2n + \)nT)], (2.43) 
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where ^ ^ 

Zb{oj) = detDo * (w)detD7'(w)detD2 - (w), (2.44) 

and 

Zf{u) = detZ)^(w)detD2(w) (2.45) 
with ^ 

Doioj) ^ --^ + co\ (2.46) 

^2 

£)i(a;) = -— + w2+e20(2 + 5), (2.47) 

^2 

D2(w) = ^ + + e^'^(4 + 7? - 25), (2.48) 

flcr2 

and ^ 

£>±(w) = — ^ +wHe^'^ + e'*4^. (2.49) 
flcr2 dcr 

Notice that D+{co) — » D-{u) upon the reflection o" — » -cr. 

The Poisson formula 

^ F(n) = ^ I t/xe^™'^"^ (x) (2.50) 
enables us to convert the infinite series (|2.43|) into [n4l 



da) 2 r°° rImlnZfc(-/nT) ImlnZK-/tir), 

rW^[C2] = - J ^ ^[lnZ,(a;) + lnZ,(c.)] + - ^ - ^J^^ \ (2.51) 

InZi(w) and lnZ/((jj) is proportional to the infinite products of (rn - E) with £"s the discrete eigenvalues 
defined by 

D^i-iE)^ = (2.52) 

for a - 0, 1,2, ±, subject to the Dirichlet boundary condition tp{+K) = 0, and these eigenvalues are all 
positive. The integral along the real axis can be carried out readily and we end up with 



Fi(r, T) - Fi(r,0) = TWlCi] - Hm ToW[C'^'] = — dojln (2.53) 

To^o ^ 171 Jo zf\u))Zf{oS) 

- ^[ - ^ Z - ^"'""^ - Z - ''''' ) - ^ Z - '''"'^ + 4 ln(l + e-P'^)] 

" £o Ei El Ey 

where E^ (A; = 0, 1, 2) stands for the common eigenvalues of Dkiu)) while £3 means the common eigenvalues 
of D±{co) and the superscript "(0)" denote the corresponding quantities at T = 0. The dimensionless function 



giirT) in (1.6) is given by 



g,(rT) = -—^[F,(r,T)-Fdr,0)] (2.54) 
with K = -1.33460... and its numerical values of the function g\{rT) will be computed in the next section. 

3. The evaluation of g\{rT) 



It follows from (2.44) and (2.45) that the quantity Z\,(pS)Zf{(jS) contains equal powers of functional determi- 
nants in the numerator and the denominator and the problem boils down to the evaluation of determinant 
ratios, which can be tackled with the Gelfand-Yaglom method |jr^[[l7|]. 



Consider two functional operators 



H,^-— + VJx), (3.1) 
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with a — 1,2, defined in the domain a < x < b under the Dirichlet boundary condition, it was shown that the 
determinant ratio 

det//2 f2(b\a) 



det//i Mb\a) 
where fa(x\a) is the solution of the homogeneous equation 



(3.2) 



Haf = 0, (3.3) 
subject to the conditions fa{0\a) - and /^(0|a) = 1. In terms of a pair of linearly independent solutions of 

T]a{a)^a(b) - T]a(b)^a(a) 



(3.3), (77„,^J 



fa{b\a) 



(3.4) 



where the Wronskian W[T]a, ^a] is x-independent. This method has been employed recently in [ p^ to calcu- 
late the one loop effective action of a Wilson line or that of a circular Wilson loop. See [ p^ for a review on 
other applications. 

Coming back to the subleading term of the heavy quark potential, the operator Ha corresponds to one 
of the operators ( 2.46 )~( 2.49). With Ha given by the operator ( 2. 46| ), eq.(3.3) can be solved analytically and 
we may choose the following pairs of independent solutions 



with their Wronskian is given by 



T]o — sinhwcr = cosh wcr. 



The equations (|33|) with //„ given by ( |2.47[ ), ( |2.48| ) and ( |2.49D 

Di(Lo)i// e2'^(2 + 6)]iff = 0, 

dcr^ 

D2(to)iff = ^ - [w^ + e^(4- + R- 26)]if/ = 0, 
do-^ 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



and 



(3.9) 



dcr^ dcr ' 

do not admit analytical solutions and we shall express the quantity of / for each case of ( 2.36 )-( 2.39 ) in 
terms of the Wronskian between a pair of special solutions. Because of the asymptotic behavior ^ p-, 
<5 - and R ^ -2 for cr close to +K, eqs.(3.7) and (3.8) have cr = +K as regular points with the same pair of 
indexes (2,-1). The equation D±(a))(p - has a regular point cr - +K with the indexes (2,-1) and cr - +K is 
an ordinary point of it. We associate rja to the vanishing solution at cr = -K and to the vanishing solution 
?Acr - K with the normalization conditions 



lim 



?/i(o') 



- lim 



?/2(o-) 



- lim 



r]+(o-) 



and 



a-^-K w2(cr -I- K)^ a^-K u^{cr + Kf <t^-K w2(cr -H 



lim 



lim 



lim 



cr'^K Cj\K - Cr)2 a^K (jj\K - Crf <r^K oj\K - af 

Furthermore, we require 

ri-{-K) = ^4K) = 



= 1 



1. 



and 



(3.10) 

(3.11) 
(3.12) 
(3.13) 
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with the prime the derivative with respect to cr. On account of the eveness of D\{(jj) and D2{co) with respect 
to (T, we have 

^i,2(cr) = 77i,2(-cr). (3.14) 
It follows from the relation between D+{ll)) and Z)_(w) that 



(3.15) 



Each differential equation of (3.7), (3.8) and (3.9) is of the form of a one dimensional Schrodinger equation 
at zero energy and does not admit a bound state subject to the Dirichlet boundary condition. Therefore we 
expect that 

Cl 2(Cl)) 



as cr — > and 



coiK - cr) 
CJ{K - cr) 

„ , , Ci,2(w) 
w(a + cr) 
C,(a;) , 

^+(cr) = — + .... 



(3.17) 



(3.18) 



(3.19) 



a){K + cr) 

as cr —> -K. The coefficients of divergence, Cxia)), C2(w) and C±{(jj) are nonzero for u + Q and are related 
to the Wronskian's via 



with j - 1,2,+. We have, in addition 



3aj 



It follows from the symmetry property ( ^.15 ) that 

C+((d) — C-{u)) = Ciiai) 



(3.20) 



(3.21) 



(3.22) 



For oj » I IK, the solutions 77's and ^'s can be approximated by WKB method and we find the asymptotic 
behaviors 



and 



where 



and 



Ci(w) = —e 



2Ka>+!^+... 



C2(U) 



3 



n ( , .\ 2Ka)+^+... 



Jo Z()S^ Jo J(\^ 



- 6sY + 3r 



^(1_^4)(,4_^4)3 



Jo ZQS^ Jo ^ s^ - ^4 



(3.23) 
(3.24) 

(3.25) 
(3.26) 



(3.27) 

The details of the derivation are shown in the appendix A. These asymptotic forms guarantees the UV 
convergence of the integral in ( 3.37 ) below. Also they will be used to check the numerical solutions. Here we 
present four curves about (3.23 ),( 3.24 ) and ( 3.25 ) at two different temperatures represented by the quantity 
s (For s - J- - tV) in Figure 1. 
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s=1.5236 



s=1 5236 




Figure 1: The top two curves represent Cie Cie and de '^'^'^ at s=1.5236, the bottom represent the same 
quantities that at s=2.0. 



The small to behavior can be obtained by introducing an another set of solutions, normalized differently 



'7i,2,+(o-) = 



'7l,2,+(0") 



^l,2,-(fr) 



and 



rj-(cr) = 



ri-ia-) 



Defining the coefficients C's by the diverging behavior 



as cr — > and 



_ Cl,2,-(W) 

ihiA^n = — + 

A - cr 



Ci2+(w) 



(3.28) 
(3.29) 

(3.30) 
(3.31) 



K + (T 

as cr —> -K, we find that Ci_2(w) ~ ur" and Ct,{(jj) ~ as w — > 0. 

Now we apply the Gelfand-Yaglom method with -K + e < cr < K - e and designate Uaito) to the 
quantity / of (||) with a - 0, 1,2,+ corresponding to the indexes of the operators ( Z.46 )-( 2.49 ), and the 
superscript "(0)" denote the coiTesponding quantities at T = 0. We find that 

t/2(w)t/2(w) 



and 



<(«) = 



Ul{oS)Ui{a))Ul{u) 
[t/<»'(a;)]ft/f>(..)[f/fV)]i 



where the superscript "(0)" denote the corresponding quantities at T = 0. It follows from (3.5) that 

sinh2(/:- e)w 



(3.32) 



(3.33) 



(3.34) 
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The symmetry (3.15) implies that U+(a>) - U-{to). In the limit e — » 



and 



f/l,2(w) 



f/±(w) 



C3(W) 



(3.35) 
(3.36) 



The counterparts at J = have the similar formula. Substituting (3.34)~(3.36) into (2.54),we have 



F,{r,T)-F,{r,Q) 



1 r 

7!" Jo 



d(i)\\\ 



C^cfcf^ sinh(">^(2/^(0)^^) 
Cf ViC2^sinh^(2/:w) 



(-^ ^ ln(l - e-^^°) - J] ln(l - e-/^^' ) - i 2 ln(l - e"'^^^) + 4 J] ln(l + e"^^^)), (3.37) 



El 



The first term of ( p7| ) depends on Ci, C2, C3 and cf ,Cf , which can be obtained by solving the 



Schrodinger like equations (3.7)~(3.9). It follows from ( 3.23 ), ( 3.24 ) and ( 3.25 ) that the integral converges. 
Notice that the temperature dependence is entirely through s, so we can use the same method to deal with 



CxXiXi and C 



(0) 



^(0) ^(0) 
1^2 >*-3 ■ 



The potentials of the Schrodinger like equations, e^^(2 + 6), e^'^(4 + R - 26) 



and e^ + e^^ are elementary functions of z but depend cr via ( p.26| ). So at the beginning, we use the forth 
order Runge-Kutta method(RK4) to the first order differential equation 



dz 
dcr 



/(./-/()) 
' 1-/0 



(3.38) 



to find the inverse of (2.26), which maps the equal step length with respect to cr e [-K+e, K] and to unequal 



step length with respect to z e [e,zo] according to (2.26), then again use RK4 with the boundary conditions 



(3.1C) and (3.13) to solve these equations. Here we present two curves about z versus cr at two different 
temperatures in Figure 2. 



Figure 2: The left corresponds to z versus cr at s= 1.5236 with K 
s=2.0 with /<: a 1.325163zo. 



1.356203zo, the right corresponds to z versus cr at 



For eqs.(3.7) and (3.8), we take advantage of the symmetry property (3.14) and evaluate the Wronskian 
by the formula 

Wi,2{cj) = -2771,2(0)77; 2(0) (3.39) 



The coefficients Ci,2(w) follow from (3.20). For eq.(3.9) with the upper sign, we develop both 77+(cr) from 
cr - -K and ^+(cr) from a - K, evaluate their Wronskian at cr = and calculate the coefficient Ciico) 



from eq.( 3.2C ). An alternative way is to iTin the solution 77+(cr) all the way to K and calculate the Wronskian 
by eq.( 3.21 ). To avoid the rapid changes of the potential function near the singularity cr - -K, we start 
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with an analytical approximation of rji^icr) and ?7+(cr) at o" = -K + p with e << p « K and then run the 
Runge-Kutta iteration for cr > -K + p. On writing x — up, we find the approximate solutions 



and 



where 



ri\.2(-K + 6)^ 3{mi(x) + c\_^2[p(x)ui{x) + q{x)vi{x)]} 



T]+(-K + (5) = 3{mi(x) + c+[/?(x)mi(x) + q(x)vi(x)]} 



ui(oJcr) — cosh WO" - 



sinh djcr 



vi(ajcr) = 1 H 



cjcr 



p(x) 



and 



20w4L3 

1 ri 3 



3 5 /I 2 3 
JC - JC H -X + -x + 



q(x) = rf-x^' - X H — ix^ + -|sinh2x- -xcosh2xl. 

^ 20w4L3 2 \ 2/ 2 J 



(3.40) 

(3.41) 

(3.42) 
(3.43) 

(3.44) 



The coefficients c\ = 4 + 4, C2 = -16 - ^ and c+ - -1 - 4- No such a precaution is necessary for 



the solution ^+(a") and the Runge-Kutta can start right at the point cr = K. The second term of (3.37) is 



dependent on £„, which refers to the eigenvalues of Dai-iE) defined by (2.52) and they are all positive. For 
Do(-iE), we find 

£o = ^ (3.45) 



with n positive integers. The eigenvalues pertaining to other D„'s satisfies the condition 

Cai-iE) = 0. 



(3.46) 



Replacing ci/ to -E^ in eqs. (3.7), (3.8) and ([3.9|), we can obtain the coefficients Cai—iE) with the same 
RK4 method and thereby its zeros. Some low-lying eigenvalues of each Dd-iE) at diff'erent temperatures 
are presented in the Table 1 . 



s 


Eu 


£12 


E21 


£22 


£31 


£32 


1.5236 


2.10070 


3.34648 


1.47031 


3.13753 


1.55800 


2.80913 


1.55 


2.10802 


3.35507 


1.49547 


3.14910 


1.56451 


2.81671 


1.60 


2.12013 


3.36926 


1.53651 


3.16823 


1.57527 


2.82923 


1.65 


2.13032 


3.38118 


1.57056 


3.18433 


1.58435 


2.83975 


1.66 


2.13216 


3.38333 


1.57666 


3.18725 


1.58599 


2.84167 


1.67 


2.13395 


3.38541 


1.58256 


3.19006 


1.58758 


2.84350 


1.68 


2.13568 


3.38743 


1.58826 


3.19279 


1.58912 


2.84528 


1.70 


2.13897 


3.39127 


1.59909 


3.19799 


1.59206 


2.84868 


1.80 


2.15270 


3.40726 


1.64378 


3.21962 


1.60430 


2.86282 


2.00 


2.17068 


3.42816 


1.70129 


3.24794 


1.62036 


2.88132 


2.50 


2.19042 


3.45103 


1.76318 


3.27901 


1.63802 


2.90159 


10.00 


2.20378 


3.46646 


1.80438 


3.29997 


1.64998 


2.91527 


200.00 


2.20445 


3.46652 


1.80621 


3.30005 


1.65214 


2.91533 



Table 1: The first and second eigenvalue of each D„(-iE) at different s(a measure of temperature). Where £„i stands 
for the first eigenvalue of D„ and £0,2 means the second eigenvalue of £)„. 

We find that the spacing between successive eigenvalues converges quickly to the asymptotic value ^ 
for large eigenvalues. Furthermore the eigenvalues hardly change when s > 2.5. Such a limiting behavior 
implies above the critical s, they almost return to the case of zero temperature. 
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In the last step, substituting ( 3.37 ) into ( 2.54 ), we finally obtain the numerical values of gi(rT) and the 
correction to the screening radius. Figure 3 is the curve of g\{rT) versus rT. The curve of goirT) is also 
displayed for reference. The physical meaning of gi{rT) will be discussed in the next section. 



1.01 -, 




Figure 3: The left curve represents gi(rT), while the right represents go(rT). 



4. Concluding remarks 

In this paper, we generalized the previous calculation of the order correction of the heavy quark 

potential at T = to the case at T > 0. We start from the partition function of the fluctuation around 
the world sheet underlying the leading order heavy quark potential with a Euclidean time Q < t < l/T 
and calculate its logarithm, the one loop effective action. The partition function is a product of ratios of 
the functional determinants of the world sheet Laplacian plus mass terms. Each determinant ratio is then 
decomposed into an infinite product over Matsubara energies with each factor a ratio of one-dimensional 
functional determinants. The summation over Matsubara energies in the eff'ective action is converted into 
an integral over Euclidean energies and an infinite series over the eigenvalues of the functional operators 
involved. The determinant ratios in the integrand are evaluated numerically by the Gelfand-Yaglom method 
and the eigenvalues are determined similarly. 

Like the zero temperature case, the order C> j correction reduces the magnitude of the heavy quark 
potential and leads to a smaller screening radius. The magnitude of reduction, measured by Kg\{rT), follows 
closely that at zero temperature within the screening radius. The function g\{rT) decreases from one at 
rT = to about 0.92 at the screening radius. Therefore, within ten percent of accuracy, the subleading 
correction may be approximated by that at T = 0, i.e. 

,,,,Ty-.^yj^ ,4.1, 

In the potential model of the heavy quarkonia, the order O j correction will lower the meson melting 
temperature. 

In the weak coupling, the heavy quark potential at T > is of Yukawa type, that is nonvanishing 
for arbitrarily large r. To the leading order of strong coupling, on the other hand, the magnitude of the 
potential drops to zero at a finite r. For a smooth cross over from the strong coupling to the weak coupling, 
one would expect the order term increases the screening radius. This would require the function 

g\{rT) vanishes somewhere within the screening radius. This, however, is not the case as shown in Figure 
4. Therefore, higher order terms are required to reverse the trend when A is lowered sufficiently. Because of 
the continuity of the free energy F{r, T) with respect to its arguments and the coupling, it is expected that 
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the cutoff nature of the potential, F(r, T) - for r > ro, remains in the weak coupling regime with ro much 
longer than the Debye length of the Yukawa potential. 




Figure 4: The left figure represents heavy quark potentials at zero temperature, while the right one represents the 
potentials at finite temperature. In both figures, the top lines correspond to the potential including subleading order 
correction and the bottom ones are leading order results. 
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A. 

In this appendix, we present the details of the derivation of Ci , C2, C3 at large a>. 
A.l Bosonic contribution 

We start with the 2nd order ordinary differential equations 

^2 

[-^+u^+e\2 + S)]m^Q (A.l) 



and 



^2 

[ ^ + + ^20(4 26)]i]2 = (A.2) 

iicr2 



As cr ^ -K, denote cr - -K + x with x «i 1, we have 

(f , 2 



(--^ + + = (A.3) 

ax x' 



2> 

(f 2 



(-^2 +^'' + = (A.4) 

ax^ x^ 



The vanishing solution at cr = -K reads 



sinhiwx^ n7 

rji = TI2 = 3(cosh(yx ) = (a>x) (A. 5) 

a>x 



-14- 



for tijx <s 1 and 

7,1 = ^2 = - — ) (A.6) 

2 ojx 



for w;ic » 1. Similarly, the approximate solution near cr = K is obtained by replacing x in (A..5) and (A.6) 
hy y = K - cr. 

Then, we have the WKB approximation for ll) '» l/K 

and 

Both WKB solutions are valid for oj{cr + K) » I and cr + K « K, and can be matched with ( |A.6| ) there. 
We have, up to an additive constant 



do-'e^^'(2 + 6') ^ I do-' e^'^' i-R' + R' + 2 + 6' 

[ dcr'-^ + [ do-'e^'(R' + 2 + (5') = 2— + f do-'(R' +2 + 6' 
J da-'^ J-K do- J_K 



d(h 2 
= 2-f = — . (A.9) 

da X 

We find that 

A = '^e'^^ (A. 10) 

Similarly 

( dcr'e^'{A + R-2S)^2— +2 \ da-' e^' {R' + 2 + 5') = 2— ^ -- (A.ll) 
J do- j^i^ do- X 



and 



3 

B = -e"'^ (A. 12) 



For K- a- <s: KbvX co{K - cr) » 1, we have 



'11 - 2^ 

^ ^^aw^+^LX^rfcre^^R+a+tJjg-dj.v^^ ^. J_-) (A.13) 

2 



where y s - cr. This is to be matched to the approximate solution 



771 oc g-'^'Xl + — ) (A.14) 



coy 



there, which is valid for o- K , i.e. y — > 0^. 
Define the coefficient Ci such that 

Ci 



We have 

for w » 1 . 
Similarly 



771-— ,(y^0+) (A. 15) 

y 



2 
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Then 



with 



772 = y,(3;^0+) 



C2 s -e^w^^+s 
for a>K » 1. 

A.2 Fermionic contribution 

The process is similar to the bosonic contribution part, the equation reads 
where 



Z)± = :r + + e^'^ + 

dcr^ 



da 



Owning to the symmetry property, we just consider 

D^ilf^ = 

For x = a + K <^ K,\he. vanishing solution at cr = -K reads 

sinhojx^ , ^2 
7]+ = 3(cosha)x ) = {(OX) 



for ojx I and 



cox 



^ le'^W - — ) 
2 ojx 



for » 1 

Similarly, the WKB approximation for a» » \IK 



Define C+ such that 
so we obtain that 

In terms of 



n^iK) = 3C3 



C3 = ^e^''-^^LU-^'''(^'^^\(aj:^ 1) 



1= f''dzyfg(R + 2) = ]- f do-e^^(R + 2) 
Jo 2 J_K 

Jrzo 1 pK 

dzy/gS=- da-e^d 
2 J_jf 



we have for w » 1 IK 



C3(w) = le2i:6.+i+..._ 
2 



(A. 18) 
(A.19) 



(A.20) 
(A.21) 

(A.22) 

(A.23) 
(A.24) 



(A.25) 

(A.26) 
(A.27) 

(A.28) 
(A.29) 

(A.30) 
(A.31) 
(A.32) 
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